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Abstract 


The  overflow  stream  from  an  M/PH/1  queue  of  finite  capacity  is  used 
as  the  input  to  an  unbounded  queue  with  one  or  more  exponential  servers. 

It  is  shown  that  the  combined  system,  consisting  of  the  two  queues,  may  be 
studied  as  a  highly  structured  Markov  process.  In  the  stable  case,  this 
Markov  process  has  a  matrix-geometric  invariant  vector.  Particular 
features  of  the  infinitesimal  generator  of  the  process  may  be  used  to  sim¬ 
plify  the  numerical  computation  of  various  steady-state  features  of  the 
model.  Several  variants  and  numerical  examples  are  discussed. 


Overflow  in  queues,  matrix-geometric  solution,  phase  type  distributions, 
systems  of  queues,  computational  probability. 


1.  The  Model 


We  study  a  service  system  with  two  units,  each  with  a  single  server. 

Unit  1  is  of  finite  capacity  K.  Unit  2  has  as  Its  Input  the  stream  of 
overflow  customers  from  Unit  1,  i.e.  any  customer  who  upon  arrival  to  the 
system  finds  K  customers  (waiting  or  In  service)  in  Unit  1  proceeds  for 
service  to  Unit  2.  The  queue  in  Unit  2  Is  unbounded.  Tractable  variants 
of  this  model,  involving  multi-server  units  with  exponential  servers,  are 
discussed  in  Sections  4  and  5. 

The  arrivals  to  the  system  are  according  to  a  Poisson  process  of  rate  X. 
The  successive  service  times  in  Unit  1  are  independent  and  have  a  common 


distribution  of  phase  type  [6,11],  given  by  the  irreducible  representation 
(j$,S)  of  order  m.  Service  times  are  assumed  to  be  positive  random  varia¬ 


bles,  so  that  g^|  ^ 


0.  Customers,  who  upon  arrival  find  K  customers  in 


Unit  1,  go  to  Unit  2.  Service  times  in  that  unit  are  independent  exponential 


random  variables  with  parameter  4.  The  interarrlval  times  and  the  successive 
service  times  in  both  units  form  independent  sequences  of  random  variables. 

For  the  present,  we  assume  that  there  are  no  independent  arrivals  directly 


to  Unit  2. 


This  queueing  system  may  be  viewed  as  an  M/M/1  queue  in  a  Markovian 
environment  [8,9,11].  In  order  to  see  this,  we  first  consider  Unit  1  alone. 


It  may  be  studied  as  a  continuous-parameter  Markov  chain  on  the  state  space 
■  {0}u{(i,j):  1  <  i  <  K,  1  <  ]  <_m}.  The  index  i  denotes  the  number 

of  customers  in  Unit  1  and  the  index  j  gives  the  phase  of  the  current  ser¬ 


vice.  The  state  0  corresponds  to  the  case  where  Unit  1  is  empty. 

The  Infinitesimal  generator  Q  of  this  Markov  process  is  given  by 
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(1) 


-X 

S° 


Q  - 


\& 

S-XI  XI 

S9B#  S-XI  XI 

S9B’  S-XI  XI 
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S*B°  S-XI  XI 
S°B°  S 

where  the  matrix  S°B®  is  defined  by  A"*A*  The  states  (i,j),  1  <  i  <  R, 
1  <_  j  <_  m,  are  ordered  in  the  lexicographic  manner. 

The  process  Q  now  clearly  acts  as  an  environmental  process,  which 
modulates  the  arrival  rate  to  the  second  unit.  Whenever  the  process  Q 
is  in  one  of  its  states  (K,j),  1  £  j  <_  m,  there  are  arrivals  to  Unit  2 
according  to  a  Poisson  process  of  rate  X.  At  all  other  times  the  arrival 
process  to  that  unit  is  interrupted. 

For  future  use,  we  note  the  following  explicit  formulas  for  the  unique 
stationary  probability  vector  i  of  Q,  which  satisfies  the  equations 
jr  Q  »  J3,  it  e_  •  1. 

Lemma  1 

After  partitioning  into  the  form  •  •  •  *Ir]  *  t*ie  vector  jj.  is 

given  by 


K-l 

-  1  +  £  Z  Cve 

__  v-1 


-i-l 


+  X 


A  CK"1(-S"1)e 


(2) 


ir^  -  iTq  A  c*  *  for  1  <_  i  £  K-l  » 

I*  -  X  irQ  A  CK_1(-S_1)  , 
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where 

(3)  C  -  A(AI-AB#*-S)"1  . 

The  matrix  B##  Is  defined  by  B#a  »  _e  •  ji  . 

Proof 

The  vector  jt  Is  the  stationary  probability  vector  for  the  (Km+1)- 
state  Markov  process  which  describes  the  bounded  M/PH/1  queue.  The  simple 
proof  of  the  lemma,  which  is  given  in  [11],  Is  repeated  here  for  completeness. 
The  steady-state  equations  jr  Q  »  may  be  written  as 

-A  ttq  +  ir  -  0  , 

A  *0  -  +  -1  (S~XI)  +  —2  S°B°  “  -2.  * 

(L) 

xLl.1  +  li^S-XI)  +  ^  S°B°  -  0  ,  for  2  <  1  ±  K-l 

A  S  ■  0  . 

Postmult ip lying  all  equations  in  (4),  but  the  first,  by  &,  we  obtain 
A  ir^  is  *  JS°,  for  1  <_  1  <_  K-l  ,  or  equivalently  A  ir^  B°°  • 

Ll+i  S*B°.  Replacing  the  rightmost  term  In  the  equations  for  1*1,  ...» 

K-l  by  A  B°°,  we  obtain 

x  ^0  -  "  -1  (XI“XB**"S) 

X  ^i-1  “  li  (XX-XB8°-S)  ,  for  2  <_  1  <  K-l  . 

This  readily  leads  to  the  formulas  (2).  The  non-singularity  of  the  matrix 
C  may  easily  be  proved  by  contradiction. 


5 


The  entire  system  may  now  be  studied  as  a  Markov  process  on  the  state 
space  E  *  {k  ^  0}  E^,  with  generator  Q*  given  by 

Q-XA  XA 

Hi  Q-XA-pI  XA 

Hi  Q-XA-uI  XA 

Hi  Q-XA-hI 

•  • 

•  • 

•  • 

where  A  denotes  a  diagonal  matrix  of  order  KnrH,  whose  last  m  diagonal 
elements  are  equal  to  one.  All  its  other  elements  are  zero.  In  Formula  (5) 
1  denotes  an  identity  matrix  of  order  Kmfl  . 

The  matrix  Q*  is  a  particular  case  of  the  generators,  discussed  in 
[11].  The  following  theorem  is  therefore  immediate. 

Theorem  1 

The  queueing  system  is  stable  if  and  only  if 

(6)  X  ®.  <  h  • 

The  stationary  probability  vector  x  of  Q*,  partitioned  into  vectors 
x^,  k  >_  0,  of  dimension  Km+1,  is  given  by 

(7)  “  I.  (I_R)  »  for  k  ^  0  . 

The  matrix  R  is  the  minimal  nonnegative  solution  to  the  equation 

(8)  H  R2  +  R  (Q-XA-hI)  +  XA  -  0  , 
and  has  spectral  radius  less  than  one. 
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As  was  pointed  out  in  [11] ,  the  fact  that  all  but  the  last  m  rows 
of  A  vanish  implies  that  the  same  is  true  of  the  matrix  R.  Although  the 
matrix  R  is  of  order  Kmfl,  only  its  last  m  rows  are  positive.  This 
special  structural  property  may  be  exploited  to  compute  R  efficiently. 

2.  Algorithmic  Procedure 

Taking  the  simplified  structure  of  the  matrices  R  and  Q  into  ac¬ 
count,  we  may  rewrite  Equation  (8)  as  follows.  We  partition  the  last  m 
rows  of  R  into  [Rq,  R^,  R2,  ....  R^] ,  where  Rq  is  a  column  vector  of 
dimension  m  and  R^,  1  <_  i  <_  K,  are  square  matrices  of  order  m.  It  then 
follows  from  (1)  and  (8)  that 

V  R^Rq  -  (X+p)  Rq  +  Rjl*  -  0  , 

u  +  XRqB  +  Rj_  [S-(X+v)I]  +  R2  S°B°  -  0  , 

P  RKRi  +  XR^  +  R±  [S-(X+p)I]  +  Ri+1  S°B°  -  0  , 

for  2  <_  i  <_  K-l  , 

P  R^  +  X  Rk_1  +  Rj^  [S-(X+p)I]  +  XI  -  0  . 

M  ■  [(X+u)I-S]  \  the  equations  (9)  may  be  written  In  the 
convenient  form 

Rq  -  (X-hi)-1  (p  . 

Rx  -  [y  skr1  +  (XRq+R210)^]  M  , 

(10)  Rt  -  [p  R^R^+XR1_1]  M  +  Ri+1  S°*0  M  , 

for  2  <_  i  <_  K-l  , 

RK  -  [Vi  R^+XR^_X+XI]  M  . 


(9) 


Setting 
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It  may  be  verified  by  an  elementary  induction  argument  that  the  Gauss- 
Seidel  type  Iterates,  defined  for  n  >_0,  by 


R^n+l)  -  (A+p)"1  [p  ^(n+l)  Rq  (n)  +  R^n+DS")  , 

Rj.Cn+l)  -  {p  R^n+1)  R^n)  +  [A  R^n)  +  R2(n+l)S°]  J3}M  , 

R^n+1)  -  [p  RK(n+l)  R^n)  4-  A  R^O*)]  M  +  Ri+1(rrfl)£°-e.  M  , 

for  2  <_  i  K-l  , 

R^n+1)  -  [p  R^(n)  +  A  R^fr)  +  AI]  M  , 

with  R^(0)  “AM,  R^(0)  ■  0,  1  £  i  £  K-l  ,  Rq(0)  ■  0_  ,  converges  increasingly 

to  the  desired  matrix  R.  This  method  of  computation  has  clear  advantages  over 
successive  substitutions. 

As  was  proved  in  [8],  the  matrix  R  satisfies  the  relation 


(11)  p  R  £  “  p 

which  provides  a  useful  accuracy  check  on  the  computation  of  R.  The  matrix 
R  is  formed  by  the  last  m  rows  of  R. 

From  Equation  (7),  by  using  the  structure  of  R,  we  obtain 

x*  jr  -  7T  R  , 

(12)  k_i  „ 

Xjj  -  Ir  (I-Rjj)  rJ  R  ,  for  k  £  1  . 


5o  + 


K 

I 

i-1 


Rje 


A  e 


By  virtue  of  (11) ,  the  marginal  density  of  the  queue  length  in  Unit  2  is 
given  by 


(13) 


yo  "  1  •  a  p'1  Ir  «  » 

yk  ■  A  p-1  (I-Rk)  R^_1.£  »  for  k  >_  1. 
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It  is  of  particular  interest  to  study  the  conditional  densities  of  the 
queue  length  in  Unit  2,  given  the  number  of  customers  present  in  Unit  1. 

To  this  end,  we  partition  the  vector  x^  as  [x£^  .x^1^ , . . .  ,x£K^  ] ,  where 
the  x£r^ ,  1  <  r  <  K,  are  m-vectors.  The  conditional  queue  length  density, 
given  that  there  are  j  customers  in  Unit  1  is  then  given  by 


(14) 


.00 


“0 


(k) 


e  , 


1  *(0) 
0  \ 


IT 


for  1<  J  <  K  , 


for  k  ^  0  . 

From  the  formulas  (12),  (13)  and  (14),  expressions  for  the  first  and 
second  moments  of  the  corresponding  queue  lengths  may  be  routinely  calculated. 


3.  The  Waiting  Time  Distribution 

We  now  consider  the  stationary  distribution  of  the  waiting  time  of  an 
arriving  customer,  under  the  f lrst-come,  first-served  queue  discipline.  There 
are  two  cases.  Either  the  customer  is  admitted  into  Unit  1  and  is  served 
immediately  or  he  waits  for  a  time  not  exceeding  x  units  of  time;  or  the 
customer  finds  the  first  unit  full  and  goes  to  the  second  unit.  If  it  is 
free,  he  is  immediately  taken  into  service.  If  not,  he  may  have  to  wait 
there  for  not  more  than  x  units  of  time.  We  denote  the  conditional  waiting 
time  distributions,  corresponding  to  these  two  alternatives  by  W^(*)  and 
W2(.)  respectively . 

W1 (x)  may  be  viewed  as  the  distribution  of  the  time  till  absorption  in 
the  Markov  Process  with  infinitesimal  generator 
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<3  0^ 

S 

S°B°  S 

0  S°B°  s 

0  S*B° 

i  «  » 

»  •  • 

S  0 
S  °B°  s. 

with  the  initial  probability  vector, 

(15)  Xq  *  (l-i^  e)"L  (ir0,  ir  ,  ir2,  ....  ir^_  )  . 

In  order  to  compute  this  distribution,  we  consider  the  system  of  dif¬ 
ferential  equations, 

y’(x)  -  Zl(x)  s°  , 

(16)  l^(x)  =  .^(x)  S  +  Xi+jOt) §.*•£.  »  for  i  <_  i  <_  K-2  , 

X^_i<x)  -Xr^Cx)  S  , 

with  initial  conditions,  given  by  (15)  . 

Clearly  [x^(x)]j  is  the  probability,  that  the  customer  is  in  the 

state  (i,j),  1  <_  j  <  m,  of  the  Markov  process  at  time  x.  From  this 

it  follows  that  W^(x)  *  yg(x),  for  x  >  0  . 

The  conditional  distribution  W^(’)  is  obtained  as  follows.  Since  the 

service  time  in  Unit  2  is  exponential  with  parameter  y,  the  Laplace-Stieltjes 

ft 

transform  W2(*)  is  given  by. 
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»;<,>  -  jo  (jj)r  x™  « 

(17)  -  Jo  [^)r  (I-V  s£  . 

-  (x^e)"1  ^  (I-R^e  +  (i^e)'1  u  %  (I-I^)  ^  [sI-uO^-I)]'1  e  . 

Upon  Inversion,  we  obtain 

(18)  W2  (x)  -  1  -  O^e)-1  itg  expfyxCR^-I)  ]e  ,  for  x  >  0  . 

To  compute  W2(*)«  it  suffices  to  solve  the  system  of  differential 
equations 

v*  (x)  -  u  v(x)  (I^-I)  , 

with  initial  conditions  v(0)  *  (x^e)  1  ,  and  to  form  W2(x)  * 

1  -  v(x)<e  ,  for  x  >_  0  . 

Finally,  the  unconditional  waiting  time  distribution  is  obtained  as, 

W(x)  -  (1-x^e)  W-^x)  +  e  W2(x)  ,  for  x  >_  0  . 

The  mean  waiting  time  of  an  arriving  customer  is  computed  by  considering 
the  two  units  separately.  The  Laplace-Stieltjes  transform  of  W^(x)  is 
given  by, 

*  K_1  -1  r 

W-  (s)  -  it-  +  2  x  [  (sI-S)  x  S°B°]  e 

1  u  r-1 

K-1  -1  r-1 

-  xn  +  2  x  (si-s)  x  s  °  r  X(s)  , 
u  r-1  ^ 

where  f(s)  -  (sI-S)  ^  S^°.  The  mean  of  W^(»)  is  hence  given  by  , 

,  K-l  .  K-l 

-  W*  (0)  -  -2  x  S'  e  +  u'  2  (r-1)  x  e  , 

1  r-1  ~ r  1  r-2  “r  “ 


where  is  Che  mean  service  time  in  Unit  1.  Similarly*  upon  differentiating 

W2(s)  we  get 

-  wj'(O)  «  u"1  (Clfci)'1  \  (I-Rk)'1e-lj  . 

The  mean  waiting  time  U  Is  then  given  by 

W«  (1-^e)  (-W*'(0))  +  (^e)  <-W*'(0))  . 

4.  Several  Servers  in  Unit  2 

Here  we  allow  the  second  unit  to  have  c  exponential  servers  of  service 
rate  u.  The  queueing  system  is  now  studied  as  a  Markov  process  on  the  state 
space  E  ■  {k  >_  0}  with  infinitesimal  generator 

A00  A01 

A10  An  Ai2 
•  •  • 

Ac-2,0  Ac~2,1  Ac-2,2 

Ac-1,0  Ac-l,l  A0 

A2  Al  A0 

^2  • • • 


The  square  blocks  of  dimension  Ko+1  are  given  by. 
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Aqo-Q-aa  . 

A01  ■  A12  *  * •*  "  Ac-1, 2  “  A0  ■  XA  ’ 

■  lyl  ,  for  1  <  i  <  c-1  , 

•  Q-XA-iyI  ,  for  1  <_  i  <_  c-1  , 

A^  “  Q-XA-cyl  , 

A2  -  cyl  , 

where  Che  matrices  Q  and  A  are  as  given  in  Section  1. 

Then,  corresponding  to  Theorem  1,  we  have 
Theorem  2.  The  queueing  system  is  stable  if  and  only  if 

X  £  <  c  y  . 

The  stationary  probability  vector  x  is  partitioned  as 

(x_,  x, ,  x„,  ...»  x  ,  x  ...)  ,  where  each  of  the  component  vectors  is  of 
— U  — X  — i  — c  — c-t-x 

dimension  Kurt-1.  The  vectors  x^,  for  k  :>  c  ,  are  given  by  , 

Cl9)  "  — c— 1  Rk’C+1  ‘ 

The  matrix  R  is  the  unique  solution  in  the  set  of  nonnegative  matrices 
of  spectral  radius  less  than  one,  of  the  matrix-quadratic  equation, 

(20)  c  y  R2  +  R  (Q-XA-cyl)  +  XA  -  0  . 

The  last  m  rows  of  R  are  positive;  all  others  are  zero. 

The  vector  (Xq,  j^,  ...,  x£_1)  is  the  left  eigenvector,  corresponding 
to  the  eigenvalue  zero,  of  the  irreducible  semi-stable  matrix 


fr* fl.rt  rl  .I  ni rift  ‘jnr  i  Y 1 1  •  • 
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A00  A01  0 
A10  A11  A12 


(21)  T  - 


A20  A21 


Ac-2,0  Ac-2,1  Ac-2,2 


A 


c-1,0 


Ac-l,l+RA2 


and  is  normalized  by 


41  -1 
2  x,  e  +  x  .  (X-R)  <e  -  1  . 

1-0  ~  ~C~L 


Proof.  The  statements  about  the  steady-state  vector  x  were  proved  In  [9] . 
That  the  matrix  R  has  the  particular  structure  Is  clear  from  the  discussion 
In  Section  1. 

The  matrix  R  is  computed  using  the  algorithmic  procedure  described  in 
Section  2,  as  are  the  steady-state  vector  and  the  other  quantities  of  interest. 


The.  Wait  in  g^ Time_Distributipn 

In  the  computation  of  the  waiting  time  distribution  of  an  arriving  cus¬ 
tomer,  only  the  distribution  W2(*)  requires  some  discussion.  Its  Laplace- 
Stieltjes  transform  W^ts)  is  given  by 


W2*8*  "  (ij.e.)"1 


C_1  (K) 

2  x^'e  +  Z 


i-0 


i-c 


-£H_) 

(cp+sj 


i-c+1 


(K) 

*1  e 


Since  R^*C+^‘,  for  i  ^  c,  we  obtain 


W2(s)  *  (V)~1 


c-1 


i-0 


x^e.  +  cp  £{  Rr  (sI+cpI-cpR^)”^e 


Upon  inversion  and  simplification,  one  obtains 
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W2(x)  -  1  -  (jr^e)  1  (I-Rjj)  1  exp  [-cyxCl-R^  ]e  , 

for  x  >_  0  . 

We  see  chat  Che  distribution  W2(*)  may  again  be  computed  by  the  solution 
of  a  simple  system  of  differential  equations,  similar  to  that  obtained  in 
Section  3. 


5.  Overflows  from  a  Multi-server  Unit  1  to  a  Unit  2  with  Several  Servers 

In  this  section,  we  consider  the  system  with  Unit  1  having  r  exponential 
servers  with  parameter  y^,  and  the  second  unit  with  c  exponential  servers, 
with  parameter  y2«  A  particular  case  of  this  model,  in  which  both  units 
have  single  exponential  servers  with  the  same  rate  and  in  which  the  waiting 
room  in  unit  2  is  also  of  finite  capacity,  has  been  studied  in  [2]. 

We  describe  the  unit  1  of  the  present  model  by  a  Markov  process  on  the 
state  space  *  {0,  1,  2,  . ..,  K},  with  infinitesimal  generator 


Q3- 


0 

1 

2 

• 

r 

r+1 

K 


-X 


V1  ~X'V1 


2yx  -X-2y1  X 


ryx  -X-ry1 
*1 


ry,  -X-ry,  X 


ryl  ^1 


The  stationary  probability  vector  z  of  Q3  is  given  by 
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The  queueing  model  now  readily  leads  to  a  quasi-birth-and -death  process 
on  the  state  space  E  ■  {1  >_  0>X  E^  ,  with  generator 

A00  A01 

A10  ^1  A12 
•  •  • 

>  Ac-2,0  Ac-2,1  Ac-2,2  , 

Ac-1,0  Ac-l,l  A0 

A2  A1  A0 

A2  •  •  • 


where 


Aqo  -  Q3  -  . 

*01  '  *12  "  "•  '  *c-2,2  “  *0  "  X*1  ’ 

A^q  •  iv2I  »  for  1  <  1  <  c-1  , 

-  Q3-XA1-ip2i  ,  for  1<1<  c-1  , 

A2  -  cy^I  , 

A1  -  Q3-AA1-cw2I  , 


and  A^  Is  the  diagonal  matrix  of  order  K+l 


whose  last  diagonal  element 
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is  one  and  all  others  are  zero. 

The  following  result  Is  then  immediate. 

Theorem  3.  Under  the  equilibrium  condition 

x  <  c  »2  • 

the  stationary  probability  vector  £  of  ,  in  partitioned  form,  is  given 
by, 

(22)  x,  ■  x  ,  R*  C+^  ,  for  i  >  c  . 

— i  -c-i  — 

The  matrix  R  is  the  unique  solution  in  the  set  of  non-negative  matrices 
of  spectral  radius  less  than  one,  of  the  equation 

(23)  c  y2  R2  +  R  CQ3-cu2I-Xd1)  +  X^  -  0  . 

The  vector  •  (Xq,  x^,  ....  x^.^),  i*  obtained  as  the  left  eigenvector 

corresponding  to  the  eigenvalue  zero  of  an  irreducible,  semi-stable  matrix, 

similar  to  the  matrix  of  Equation  (21) .  It  is  normalized  by  the  condition 

1  x,  e  +  x  .  (I-R)  e.  *  1  . 
i-0  ^ 

Only  the  last  row  of  R  differs  from  zero  and  is  strictly  positive. 


The  elements  Rq,  R^»  ...»  R^  in  the  last  row  of  R  are  computed  by 
iterative  solution  of  the  system  of  equations 

cp 2  \  *0  "  <X+CW2^  R0  +  W1  R1  *  0  » 

cy 2  Rg  R j  +  X  Rj_2  -  (X+jM^+CMj)  Rj  +  (j+1)  Mj  *  0  , 

for  1<  j  <  r-1  , 
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cp2  R^Rj  +  X  -  (X+rv1+cw2)  R^  +  ryx  R^+1  -  0  , 

for  r  <  j  <  K-l  , 

cy2  R^  +  X  Rr_1  -  (X+ru1)  Rj^  +  X  -  0  . 

The  accuracy  check 

should  be  satisfied. 

6.  Additional  Variants 

For  each  of  the  three  preceding  models,  the  Unit  2  may  be  assumed  to 
have  Its  own  Poisson  arrival  process  of  rate  X  .  The  general  results  on  the 
matrix-geometric  form  of  the  stationary  probability  vector  x  remain  valid 
for  this  case,  but  the  major  simplifications  due  to  the  particular  structure 
of  R  are  lost  for  X^  >  0.  The  matrix  R  is  then  strictly  positive.  In 
the  computation  of  R,  It  Is  still  worthwhile  to  write  that  matrix  In  partitioned 
form,  so  that  the  highly  sparse  and  regular  structure  of  the  coefficient 
matrices  may  be  exploited.  How  this  may  be  done  is  obvious  once  the  analogues 
of  the  equations  (9)  are  written  down. 

The  overflow  streams  from  the  bounded  M/FH/1  and  M/M/r  queues  are  parti¬ 
cular  cases  of  the  versatile  Markovian  point  process,  introduced  in  [10].  By 
Implementing  the  general,  but  much  more  belabored  algorithm  in  [12],  it  is 
possible  to  study  the  case  where  Unit  2  has  a  single  server  with  a  general 
service  time  distribution.  If  one  allows  the  service  time  distribution  in 
that  unit  to  be  of  phase  type,  the  methods  of  the  present  paper  may  be  im¬ 
plemented.  The  large  matrices,  which  now  arise,  need  to  be  handled  with  care 
to  avoid  problems  of  storage  and  large  computation  time.  A  brief  discussion 
of  this  case  may  be  found  in  Chapter  6  of  [11]. 


K 

Z  Rj  ' 
j-0  3 


cu. 
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7.  Numerical  Examples 

In  some  of  Che  small  number  of  papers,  which  give  numerical  results  for 
overflow  models,  it  has  been  noted  that  such  results  may  be  quite  astonishing 
and  that  naive  approximations  to  the  overflow  stream  may  lead  to  gross  errors 
[3-5].  This  observation  may  be  made  generally  for  queues  in  which  arrival 
and  service  rates  are  subject  to  (random)  fluctuations  [8,  11].  The  under¬ 
lying  cause  for  the  highly  overdispersed  queue  length  and  waiting  time  dis¬ 
tributions  for  such  models  lies  in  the  large  variations,  which  typically  occur 
in  the  Interarrival  times  to  Unit  2.  In  order  to  illustrate  this  adequately, 
it  is  necessary  to  consider  several  related  numerical  examples. 

We  consider  the  system  consisting  of  an  M/PH/1  queue  of  finite  capacity 

K  ■  4,  with  overflow  to  a  single  exponential  server.  The  service  time  dis- 

38  1 

tribution  in  Unit  1  was  chosen  to  be  (8»x)  +  39'  E^CO.l.x).  Its  mean 

is  0.5.  Its  particular  form  was  chosen  to  reflect  a  case  in  which  occasionally 
long  service  times  occur.  The  arrival  rate  X  to  the  system  is  20,  so  that 
Unit  1  is  saturated  90X  of  the  time. 

The  service  rate  y  in  Unit  2  was  chosen  so  as  to  obtain  three  different 
cases,  respectively  with  p  *  0.92,  0.94  and  0.96.  Before  giving  the  numerical 
results,  that  are  obtained,  we  briefly  consider  an  appealing,  but  naive 
"approximation".  It  may  be  argued  that  the  overflow  stream  from  Unit  1  could  be 
"approximated"  by  a  Poisson  process  of  rate  0.9X.  This  would  lead,  by 
elementary  formulas  for  the  M/M/1  queue,  the  following  values  for  the  mean  and 
the  standard  deviation  of  the  queue  length  in  Unit  2. 


£ 

mean 

st. 

dev 

0.92 

11.50 

11. 

.99 

0.94 

15.67 

16. 

.16 

0.96 

24.00 

24. 

,49 
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The  correct  computed  values  are  however  much  larger. 


jO 

mean 

st.  dev. 

0.92 

18.99 

23.39 

0.94 

27.64 

33.44 

0.96 

45.24 

53.33 

The  reasons  why  the  naive  approximation  fares  so  poorly  are  similar  to  those 
discussed  in  [8].  During  periods  of  overflow,  the  second  unit  behaves  like 
a  mildly  unstable  M/M/1  queue,  which  may  recover  from  its  build-up  during 
the  periods  when  Unit  1  is  not  filled  to  capacity.  During  the  rare  long 
service  times  in  Unit  1,  the  latter  will  almost  certainly  become  saturated 
and  remain  so  for  a  fairly  long  time.  The  build-ups  in  the  second  unit  will 
then  be  substantial  and  cause  large  queue  lengths  for  a  long  time  thereafter. 

The  numerical  results  for  Unit  2  are  not  significantly  affected  by  a 
moderate  increase  in  the  capacity  of  Unit  1.  A.  change  of  e.g.  K  -  4,  to 
K  ■  8,  resulted  in  only  a  minute  change  In  the  computed  characteristics  of 
Unit  2.  The  additional  four  waiting  spaces  are  occupied  most  of  the  time 
and  do  not  appreciably  affect  the  rate  of  overflow. 

As  is  to  be  expected,  the  queue  in  Unit  2  is  much  smoother  when  the 
service  times  in  Unit  1  exhibit  less  random  variation.  This  is  seen  by 
replacing  the  earlier  service  time  distribution,  e.g.  by  E^(3,x). 

A  variety  of  other  aspects  of  overflow  in  queues  may  easily  be  numeri¬ 
cally  Investigated  by  implementing  the  algorithm  proposed  In  this  paper.  The 
numerical  results  are  often  astonishing  at  first,  but  are  seen  -  after  same 
reflection  -  to  correspond  to  intuitive  qualitative  behavior  of  the  queue . 
They  suggest  that  great  care  is  needed  in  the  Interpretation  of  numerical 
results  for  networks  with  capacity  constraints. 
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